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(57) Abstract 

This invention describes a method for aligning radi- 
ation therapy beams with a treatment target of a patient. 
Diagnostic computer tomography scanning is used to map 
the treatment target in the patient's anatomy and to aid 
physicians in deciding where to aim the radiation therapy 
beams. Digitally reconstructed radiographs are generated 
from an intermediate 3-D image produced from the CT 
scans. These DRRs are compared with x-ray images taken 
of the treatment target's position taken at the time of treat- 
ment A transformation equation is used to compare the 
DRRs and the x-ray images. A minimization equation is 
used to determine the similarity between the two sets of 
images. If the difference between the two sets fall below a 
determined minimum value, then the position of the treat- 
ment target is deemed correct and the radiation therapy 
begins. If the difference does not fall below the deter- 
mined minimum value, then either the radiation beams or 
the treatment target is adjusted and the process is repeated. 
This procedure may be used many times in the course of 
a single treatment session, depending on the movement of 
the patient. 
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High-Speed Inter-Modality Image Registration Via 
Iterative Feature Matching 



RELATED APPLICATION INFORMATION 

This application claims priority from U.S. Provisional 
Application No. 60/021,588 entitled "HIGH-SPEED INTER- 
MODALITY IMAGE REGISTRATION VIA ITERATIVE FEATURE MATCHING" 
filed July 11, 1996, which is herein incorporated by 
reference . 

FIELD OF THE INVENTION 

This invention relates to the field of medical imaging. 
More particularly, it relates to a real-time method of 
positioning therapeutic radiation beams with respect to a 
target area within a patient. 

25 BACKGROUND OF THE INVENTION 

Radiation therapy is often used to treat: cancerous 
tumors within a patient's body. An early diagnostic 
session is conducted where the physician uses an imaging 
30 technique, such as computed tomography (CT) scanning or 

magnetic resonance imaging (MRI) to study the target area. 
He or she then decides the ideal placement and volume of 
the radiation beam{s) with respect to the target area. 
During the actual treatment, the radiation beams are 
35 focused directly at the target area, using the diagnostic 
studies as a position reference. Precise positioning of 
the radiation beams insures that most of the radiation 
contacts the target cells, while also insuring that the 
healthy cells surrounding the target cells are not 
affected. Unfortunately, it is often difficult to be 
certain that radiation beams are optimally positioned with 
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respect to target cells. Often, a smaller total dose of 
radiation must be used in order to reduce the possibility 
of damage to healthy cells. The consequence, however, is 
that the radiation treatment becomes less effective. 

In addition, radiotherapy often requires a patient to 
return for treatment over the course of several days. 
Repositioning a patient precisely each time can be time- 
consuming and frustrating. 

Over the past decade, many methods have been devised 
to improve the alignment of radiation beams with the target 
area of a patient. An early method involves a rigid frame 
to physically hold in place the part of the patient's body 
to be treated. In one embodiment for treatment of target 
areas within a patient's skull, the frame is attached to a 
15 floorstand mounted in a Linac (linear accelerator) floor 
turret. This method is considered generally reliable and 
accurate, as it fixes the target: area rather precisely with 
respect to the radiation beams. Unfortunately, due to the 
nature of the frame, it also greatly limits accessibility 
20 to the patient's skull. Target areas may be located in the 
skull where the Linac radiation beams cannot reach. In 
addition, it is extremely uncomfortable for the patient, 
who must remain in an awkward position for long periods of 
time. 

Another method involves invasive techniques. Patent 
No. 5,097,839 by Allen describes fiducial implants attached 
in a pattern to a patient's skull bones, underneath the 
skin. These implants are then used as fixed references 
when aligning radiation beams with the target area. These 
implants are an improvement over rigid frames in that they 
allow all target areas within a skull to be reached. 
However, because inserting the fiducial implants into the 
patient is a surgical procedure itself, the patient must 
often wait for several days until the radiation treatment. 
35 During this time, the target area may grow or otherwise 

change in shape, rendering inaccurate the early diagnostic 
analyses taken when the fiducial implants were put in 
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place. In addition, the implants are often disfiguring and 
painful to the patient. 

Another type of invasive technique involves placing 
tattoos on the patient's skin where the radiation beams are 
5 to enter. Although this is less intrusive than the 

fiducial implants, it has many of the same problems, such 
as the patient having to wait several days from the time of 
the tattoo procedure until the radiation treatment, thus 
giving time for the target area to grow or change shape. 

10 In addition, given the nature of tattoos, it is possible 
they may also change shape. 

More recently, non- invasive, non-disfiguring alignment 
systems have been developed. These typically use signal 
processing to convert the CT or MRI data of the position of 

15 the patient in the diagnostic setting to the position of 
the patient in the treatment setting. Many of these 
systems require a large amount of preprocessing, whereby 
data generated from the diagnostic scan is gathered and 
manipulated until it is usable in the treatment setting. 

20 The preprocessing step can take several days. During 
treatment, real time images are compared with the 
preprocessing data and the patient or the radiation therapy 
beams are adjusted accordingly. Oftentimes, manual 
adjustment is necessary. Three degrees of freedom, 

2 5 corresponding to one plane, are typically allowed. The 

patient has greater freedom of movement than in the 
previously described techniques, but his movement is still 
confined. These systems are generally accurate, and 
painless for the patient. 

3 0 Patent No. 5,295,200 by Boyer et al . describes a 

method of aligning radiation therapy beams using the Fast 
Fourier Transform (FFT) to compare the position of 
diagnostic images with the position of treatment images. 
In this invention, a large amount of complex data must be 
35 gathered and processed prior to treatment. Reference 

images collected during the diagnostic study are later used 
to position the patient during treatment. 
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Patent No, 5,531,520 by Grimson et al . describes a 
method of image registration that takes into consideration 
patient movement within six degrees of freedom. It employs 
lasers to determine patient position and as such is 
5 confined to surface images. Thus, treatment beams must be 
based relative to tattoos or other markers on a patient's 
skin, which have the problems mentioned above. 

These existing alignment methods require an extensive 
amount of time to process complex diagnostic data, usually 

io restrict accuracy to three degrees of freedom, limit 
patient movement, and make adjustment of either the 
treatment beams or the patient difficult. In addition, 
they are unable to generate instant reference images with 
which to compare the present position of a patient. They 

15 also require manual operations to supplement automatic 
procedures . 

OBJECTS AND ADVANTAGES OF THE INVENTION 

20 Accordingly, it is a primary object of the present 

invention to provide a wholly automatic method of aligning 
radiation therapy beams with the treatment target of a 
patient. It is another object of the invention to provide 
such a method by making use of radiographic techniques . 

25 Yet another object of the invention is to decrease the time 
required for the preprocessing step, as well as reduce the 
complexity of data manipulation. A further object of this 
invention is to generate instant reference images derived 
from the diagnostic study with which to compare the present 

30 position of the treatment target. Another object of the 
invention is to allow a patient to move around freely 
during treatment. A sixth object of the invention is to 
measure patient movement within six degrees of freedom. A 
seventh object of the invention is to provide continuous 

35 adjustment of radiation therapy beams to improve precision 
during the course of treatment. A final object of this 
invention is to accommodate both diagnostic imaging and 
treatment of a patient in the same day. 
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SUMMARY OF THE INVENTION 

These objects and advantages are attained by the 
present invention. This method begins with a diagnostic 
5 computed tomographic (CT) scan of the treatment target of a 
patient. The information from the CT scan is used to 
generate an intermediate 3-D image. This intermediate 3-D 
image will eventually be moved or rotated to mimic the 
position of the treatment target. During treatment, the 

10 position of the treatment target in relation to the 

radiation therapy beams is recorded using at least two x- 
ray cameras. Thus, both translational and rotational 
information is received. The treatment images are then 
processed to produce a feature vector, which is specific 

15 for those set of images and thus the position of the 
treatment target. 

Using the intermediate 3-D image generated during the 
diagnostic stage, at least two digitally reconstructed 
radiographs (DRRs) are produced. These radiographs are 

20 wholly artificial constructs that are judged to be similar 
to the treatment images. These DRRs are also processed to 
produce a feature vector. The difference between the two 
feature vectors is then calculated using a mathematical 
equation, for example chi squared. If the difference falls 

25 below a minimum allowable value, then the position of the 
treatment target is deemed accurate enough to begin 
radiation therapy. However, if the difference does not 
fall below the minimum allowable value, then the 
intermediate 3-D image is iteratively moved within six 

30 degrees of freedom until successive DRRs produce an 
acceptable minimum difference with respect to the treatment 
images. Data from the repositioned intermediate 3-D image 
is used to adjust either the patient or the radiation 
therapy beams, thus achieving accurate alignment of the 

35 radiation beams and the treatment target. 

DESCRIPTION OF THE FIGURES 
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Fig. 1 is a flow chart illustrating the operation of the 
invention . 

Fig. 2 shows a radiograph produced from a patient with a 
brain tumor. 

5 Fig. 3 is a diagram of the apparatus, consisting of a 
Cyberknife connected to a computer system. 
Fig. 4 is a diagram showing the diagnostic coordinate 
system. 

Fig. 5 is a diagram showing the treatment coordinate 
10 system. 

Fig. 6 is a diagram showing the three translational degrees 

of freedom in the preferred embodiment. 
Fig. 7 is a diagram showing the three rotational degrees of 

freedom in the preferred embodiment. 
15 Fig. 8 is a diagram illustrating the projection geometry 

for generating two digitally reconstructed radiographs 

(DRRs) from the intermediate 3-D image using a 

hypothetical camera model. 
Fig. 9 shows a fluoroscopic image of an anthropomorphic 
20 phantom, a DRR of said anthropomorphic phantom, and a 

DRR marked with key pixels to locate anatomical edge 

features . 

Fig. 10 contains graphs showing the empirical results for 
determining (x, y, z) translations, (a,p\Y) rotations 
one axis at a time, and composite (ot,p) rotations for 
an anthropomorphic phantom. 
Fig. 11 is a table showing the standard deviation of errors 
in measuring translations and rotations of the 
anthropomorphic phantom. 
Fig. 12 is a graph illustrating the distribution of 

rotation measurement errors for the empirical tests. 
The distribution is approximately Gaussian with zero 
mean. 

Fig. 13 is a graph showing the distribution of %2 from the 
35 empirical tests. 

Fig. 14 is a graph showing the correlation between %2 and 

the errors in angle measurements from the empirical 
tests . 
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Fig. 15 is a graph showing the correlation between %2 and 
the rotation angles around each axis from the 
empirical tests. 

5 

DETAILED DESCRIPTION 

For illustrative purposes, the method of the invention 
is described for use with the Cyberknife, as shown in Fig. 

10 3, and a patient 30 with a brain tumor 32 , as shown in 
Fig. 2. Please note the method of this invention is not 
confined for use with the Cyberknife. In the preprocessing 
step, patient 30 undergoes a computed tomography (CT) scan 
of his skull. The CT scans are used to assemble an 

15 intermediate 3-D image 54, as shown in Fig. 8. 
Intermediate 3-D image 54 can be moved and rotated along 
both trans lational axes 50 and rotational axes 52 in the 
diagnostic coordinate system, as shown in Figs. 6 and 7. 
In this representation, translat ional axes 50 are 

20 represented by dx, dy, and dz, while the rotational axes 52 
are represented by a, P, y. Thus up to six degrees of 

freedom of movement are allowed. Intermediate 3-D image 54 
may consist solely of computer script within a computer or 
system of computers 44 or it may be visualized on a 
25 monitor. 

In the preferred embodiment, intermediate 3-D image 
54 is used to generate a set of at least two digitally 
reconstructed radiographs { DRRs ) 56, as shown in Fig. 8. 
DRRs 56 are artificial 2-D images that show how 

30 intermediate 3-D image 54 would appear from different 
angles using a hypothetical camera model 58. Each DRR 56 
of a set represents one hypothetical camera angle. These 
DRRs 56 can then be masked to isolate key pixels associated 
with anatomical edge features 64, as shown in Fig. 9. 

35 Where a complete image 60 would have 40,000 pixels, for 
example, a masked image 62 typically has 1,000 to 4,000 
key pixels 64 . 

The set of masked DRRs 62 may then be used to 
generate a lookup table. The lookup table provides the 



7 



WO 98/02091 



PCT/US97/11592 



first derivatives of said translational and rotational 
measurements of intermediate 3-D image 54. These 
calculations can be used later during treatment to match 
the actual patient images with intermediate 3-D image 54. 
The preprocessing procedure as described above requires 
about 20 seconds on a computer or computer systems 44 with 
200 MHz processors. 

During treatment, the patient is placed within the 
view of at least two radiographic cameras 42, in a position 
approximating his position during the diagnostic scan. The 
patient has complete freedom of movement, as all possible 
positions can be defined within six degrees of freedom 
using translational axes 50 and rotational axes 52, as 
shown in Figs. 6 and 7. Translational axes 50 and 
rotational axes 52 of the treatment coordinate system, as 
shown in Fig. 5, are defined in the same manner as 
translational axes 50 and rotational axes 52 of the 
diagnostic coordinate system, as shown in Fig. 4. (On the 
Cyberknife, patient 30 has two x-ray cameras 42 and 
screens 36, which produce real-time radiographs 31 of 
treatment target 32). These real-time radiographs 31 may 
then be processed in the same manner as DRRs 56. Real-time 
radiographs 31 are masked to isolate key pixels associated 
with anatomical edge features 64. Masked real-time 
radiographs 31 are used to produce a first feature vector, 
which specifically identifies the position and orientation 
of treatment target 32 within the treatment coordinate 
system, as shown in Fig. 5, at the time real-time 
radiographs 31 were taken. 

Next, intermediate 3-D image 54 is manipulated until 
its position emulates the position and orientation of 
treatment target 32. New DRRs 56 are then generated from 
intermediate 3-D image 54, as shown in Fig. 8. These DRRs 
56 are masked to isolate the same key pixels 64 as in 
real-time radiographs 31 and processed to produce a second 
feature vector, which specifically identifies the position 
and orientation of the treatment target of intermediate 3-D 
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image 54 within the diagnostic coordinate system, as shown 
in Fig. 4 . 

The two feature vectors are then compared using a 
mathematical equation, for example the chi squared 
5 statistic. If treatment target 32 of patient 30 is 
positioned in the treatment coordinate system, as shown in 
Fig. 5, in precisely the same way as intermediate 3-D image 
54 is positioned in the diagnostic coordinate system, as 
shown in Fig. 4, then the difference between the two 
10 feature vectors, or % 2 * will be less than a designated 

minimum value. The system has then completed its 
determination of treatment target 32 position and 
orientation. This information is passed on to the beam 
delivery system 40 (e.g. the Cyberknife) and the radiation 
is therapy beams 38 are allowed to operate. 

If, however, treatment target 32 is not positioned in 
the same position and orientation as in the diagnostic 
coordinate system, the two feature vectors will exhibit y 2 

greater than the designated minimum value. In this case, 
20 the system moves to the % 2 minimization step. The y 2 

minimization process searches for a match between real-time 
radiographs 31 and DRRs 56 by iteratively varying the 
position and orientation of intermediate 3-D image 54 with 
respect to hypothetical camera model 56. For a 

25 radiographic imaging system consisting of at least two CCD 
f luoroscopes , hypothetical camera model 58 is a multi- 
parameter function which maps the CCD image plane to the 
fluoroscope screen, and the fluoroscope screen to the 
treatment coordinate system, as shown in Fig. 5. It 

30 accounts for magnification, relative positioning of the two 
fluoroscopes in the treatment room coordinate system, tilt 
of the image planes with respect to the fluoroscope screens 
(in three independent directions), and radial (spherical) 
distortion of the lens system. Radial distortion is 

35 modeled by the factor: 

x 1 = x / (1 + XR 2 ) 
y = y / (l + Xr 2 ) 
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where x and y are the pixel coordinates on the CCD image 
plane, x' and y are the corresponding coordinates on the 
fluoroscope screen, r2 = ( X 2 + y 2 ) t and x is a free 

parameter that is determined when hypothetical camera model 
58 is calibrated. 

The parameters of hypothetical camera model 58 are 
determined by fitting DRRs 56 of a calibration phantom to 
actual radiographs 31 of the phantom, using the same % 2 
minimization process as for patient position measurements. 
The residual errors in this calibration-fitting process are 
spatially random and on the order of the image pixel 
dimensions, indicating that there is no significant 
distortion left unaccounced for by the model. 

The various positions of the treatment target are 
emulated in DRRs 56 by changing the point of view of 
hypothetical camera model 58 with respect to intermediate 
3-D image 54. This movement is achieved by use of the 
lookup table created in the preprocessing step, as 
20 described above. The process continues until the 
difference between the two feature vectors falls below the 
designated minimum value. At this point, the new 
coordinates of hypothetical camera model 5 8 in the 
diagnostic coordinate system are translated into 
25 coordinates for the treatment coordinate system. This 
information is sent to beam delivery system 40 and 
radiation therapy beams 38 are realigned. It is also 
possible to reposition patient 30. Radiation therapy beams 
38 are then allowed to operate. 
30 Translation of intermediate 3-D image 54 rotation 

geometry into the treatment target geometry is accomplished 
most effectively via an Eulerian (body-fixed) rotation 
convention. The Eulerian convention properly represents 
the fact that patient's rotational axes 52 are fixed in 
35 anatomy and not in an external reference frame. On the 
other hand, mechanical beam and patient alignment systems 
can measure angles in the space-fixed (a, p,y) convention 
illustrated in Fig. 10. To relate space-fixed rotations in 
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(a, p, y) with the rotations deduced by the image 
registration algorithm requires that the Eulerian DRR 
rotations be inverted to correspond to rotations of the 
object rather than hypothetical camera model 58 (bearing 
in mind that sequential rotations do not commute) and then 
transformed to the space- fixed rotation convention. This 
transformation is summarized below. 

Let us designate treatment target 32 in the skull of 
patient 3 0 by the vector X in the space- fixed camera 
coordinate system. When treatment target 32 rotates 
through the angles (a, p, y) , its coordinates in the fixed 

frame become 

X- - [A Y ]fAp][Aa]X 

where the rotation matrices are: 

[Ay] =1 0 0 

0 cosa sina 

0 -sina cosa 



[Apj = cosP 0 sinP 

0 10 
-sinP 0 cosp 



[A a ] = cosy -siny 0 
siny cosy 0 
0 0 1 



This convention locates points in the anatomy of 
patient 30 in the fixed treatment coordinate frame that 
beam delivery system 40 refers to in directing treatment 
beam 38. 

The DRR ray-tracing process works most effectively in 
a coordinate system fixed to the image plane in 
hypothetical camera model 58. The image plane of 
hypothetical camera model 58 is positioned within the 
diagnostic coordinate system, relative to intermediate 3-D 
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image 54, through the three Eulerian rotations ( <K), 6(), yo ) . 
Rotations of the patient's anatomy with respect to x-ray 
cameras 42 can then be represented in the DRR process by 
moving the image plane of hypothetical camera model 56 
through Q 0 , y 0 ) to the perturbed orientation «j>o + d<)>, On + 
de,yn + dy). if a point in the patient's anatomy is located 
by the vector X' in the coordinate frame fixed to the image 
plane, then X and X" are related according co the 
transformation : 

X" = [A<j>] [Ae] [Ay]X 

where the rotation matrices are: 

15 [A<j,] =1 o o 

0 cos0 sin<j> 



10 



25 



35 



0 -sin<J) cos<(> 



[Aei = cose o 



sine 



20 0 1 



0 



-sine 0 cose 



[Ay] = cos\|/ -siny 0 
sin\|/ cosy 0 
0 0 



Let [En] = [A<j>] [Ae] [Ay] define the complete Eulerian 
rotation matrix for angles (4>0. 60, VO) and [E 0 + dE] define 
the complete rotation matrix for rotations «|>o + d<j>, Q 0 + d6, yo 
30 + dy). The inverse transf ormt ion is [E~l] = 
[A^] [A_e] [A_y] . Following this convention, the vector 
locating the point X" in the treatment coordinate frame 
after a rigid-body rotation (d*. d6. dy) that is modeled by 
the rotation of hypothetical camera model 58 is given by: 



X" = [BO -1 ] [Eo + dE]X. 
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From this, the relationship between the space- fixed 
rotations ( a, P, Y) and the body-fixed Eulerian rotations 
(d<t>, d6, d\y) is given by: 

[Eo* 1 ][E 0 + dB] = [Ay][Af}][Acx] 
EMPIRICAL TESTS 

The reduction of this method to practice has been 
demonstrated in the following tests. In the tests, the 
precision of measuring actual phantom rotations will be 
reported in (a, P, y) , while numerical simulations of 
rotation measurement will be reported in (0, 8, y) Eulerian 
angles . 

In the tests, an anthropomorphic skull phantom was set 
up on the treatment couch, with the inferior/superior axis 
along the x-ray camera x-axis, as it would be for a typical 
patient. The phantom's anterior direction was in the 
positive 2 direction. The three axes of rotation in the 
test setup corresponded to the space-fixed angles (a,p, y) , 
as defined in Fig. 7. A tiltmeter measured angles (a, P) 
relative to gravity and was calibrated on a milling 
machine's rotary table. The resulting calibration was 
accurate to better than 0.1 degrees. The remaining 
translational and rotational degrees of freedom were 
measured mechanically and had precisions of 0.25 mm and 
0.25 degrees, respectively. 

The phantom was imaged in a CT study of 74 slices 3.0 
mm thick, with each image having 512 x 512 pixels 0.49 mm 
on a side. The treatment room imaging system operated at 
100 kV and 25 mA. The radiographs were 200 x 200 pixels, 
with each pixel 1.3 0 mm on a side. 

The phantom was moved step-by-step through motions 
along and around each axis separately and in composite 
multi-axis rotations. At each step, the imaging system 
acquired radiographs and the algorithm computed the phantom 
position and orientation relative to the CT study. The 
correspondence between the measured and actual change in 
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position was recorded for each individual degree of 
freedom. For each test position, the minimum value of x 2 
at convergence was recorded. 

Fig. 10 illustrates the results for determining 
translations and rotations empirically. This figure 
displays the known translations and rotations along the 
abscissa and the deduced translations and rotations along 
the ordinate. The figure shows data for rotations (a, fry) 
around one axis at a time, and data for composite rotations 
in which a and 0 were varied simultaneously. The 
individual p and y rotations each correspond to a composite 
rotation of <p and y in the Eulerian system used to model 
phantom position in the DRRs . The composite rotations 
involved all three Eulerian angles. 
15 The deduced versus actual translations and rotations 

should follow straight lines of slope equal to one. The 
data displayed in the figures were fit to straight lines, 
in one case fixing the slope to one, and in the other case 
allowing the slope to be a free parameter. The variance of 
the individual points about the fitted line is the standard 
deviation. The results of analyzing the data via straight 
line fits are summarized in Fig. 11. 

Fig. 12 plots the frequency distribution of angle 
measurement errors. This distribution has an approximately 
Gaussian shape with zero mean, which is consistent with the 
conclusion that the process is not making systematic errors 
in determining the phantom position. The frequency 
distribution for x 2 at convergence is plotted in Fig. 13. 
This distribution has the basic functional form of a X 2 
distribution for fits to data with randomly distributed 
statistical fluctuations. This again supports the 
conclusion that the process of extracting the image moments 
and edge coordinates that are fit by X 2 minimization is not 
making systematically biased determinations of the image 
35 features. 

Fig. 14 illustrates the correspondence between the 
measurement error for each rotational component and the 
value of x 2 at convergence The relationship is 
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uncorrelated for values of % 2 less than one and rotation 
errors less than about one degree. This is consistent with 
the supposition that once the position determination has 
gotten close, the minimization routine encounters a % 2 
5 surface without a sharply defined minimum. Values of j} 
greater than one have a positive correlation with 
increasing error in the angle determinations. This is a 
valuable property, as it allows one to use the magnitude of 
X 2 to flag those rare instances where a poor determination 

10 of the patient orientation has been made. 

The correlation between x 2 anc * the magnitude of 

rotation around each axis is displayed in Fig. 15. For the 
a and P rotations there is no apparent correlation, which 

indicates that the algorithm's precision and reliability 
15 is uniform over the full range of allowable orientations 
around these axes. The plot for the y rotations shows 

greater difficulty in establishing orientation in one of 
the two directions . 

20 APPLICATIONS 

The therapy beam alignment method described in the 
present invention can be used to direct radiation beams to 
any part of the body containing radiographic features. In 

25 addition, it is obvious that this method may also be used 
to align instruments other than radiation beams with 
objects other than disease regions of a patient's body. 
For example, this method would allow precise positioning of 
fabrication tools with respect to a manufactured object. 

30 Thus, it is obvious that modifications and variations 

of the present invention are possible. Therefore, it is to 
be understood that the scope of the invention should be 
based on the following claims. 

35 
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CLAIMS 

What is claimed is: 



1 1. A method of aligning a plurality of radiation 

2 therapy beams with a treatment target of a patient, 

3 comprising the steps of: 

4 a) preprocessing, comprising the steps of: 

5 1) conducting a computed tomography (CT) scan; and 

6 2) assembling an intermediate 3-D image from said 

7 scan ; and 

8 b) treatment, comprising the steps of: 

9 1) acquiring a set of at least two real images of 
10 position of said treatment target of a patient; 
n 2) using said at least two real images to produce 

12 a first feature vector; 

13 3) generating a set of at least two digitally 

1 4 reconstructed images (DRRs) from said 

15 intermediate 3-D image; 

16 4) using said set of at least two digitally 

1 7 reconstructed images to produce a second 

18 feature vector; 

!9 5) calculating a difference between said first 

20 feature vector and said second feature vector; 

21 and 

22 6) if said difference does not fall below a 

23 minimum allowable value, employing an algorithm 

24 -c adjust positions of said radiation therapy 

25 beams such that said radiation therapy beams 

26 are aligned with said treatment target of said 

27 patient. 
28 

1 2. The method of claim 1 comprising repeating steps 

2 b) 1-6. 

3 

1 3. The method of claim 1, including the step of 

2 generating key pixels associated with anatomical 

3 edge features. 
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1 4. The method of claim 3, including the step of 

2 generating a mask for said at lease two real 

3 images and said at least two digitally 

4 reconstructed images by highlighting key 

5 pixels . 
6 

1 5. The method of claim 3, wherein said feature 

2 vectors are determined by said key pixels . 

3 

1 6. The method of claim 1, using a preprocessing 

2 coordinate system and a treatment coordinate 

3 system. 

4 

1 7. The method of claim 6, using said preprocessing 

2 coordinate system wherein said treatment target 

3 of said patient is fixed. 

4 

1 8. The method of claim 6, using said treatment 

2 coordinate system wherein said treatment target 

3 of said patient is rotated and translated. 

4 

1 9. The method of claim 6 wherein said 

2 preprocessing coordinate system and said 

3 treatment coordinate system are related by a 

4 transformation equation. 

5 

1 10. The method of claim 1, allowing six degrees of 

2 freedom of position of said treatment target of 

3 said patient. 
4 

1 11. The method of claim 10, wherein said six 

2 degrees of freedom include translation 

3 measurements . 

4 

1 12. The method of claim 10, wherein said six 

2 degrees of freedom include rotational 

3 measurements , 

4 
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1 13. Trie method of claim 1, wherein said difference 

2 between said vectors is determined by use of chi 

3 squared . 

4 

1 14. The method of claim 1, wherein said algorithm 

2 provides minimization of chi squared until 

3 convergence with said minimum allowable value. 
4 

1 15. The method of claim 14, including the step of 

2 adjusting patient anatomy in said intermediate 

3 3-D image. 

4 

1 16. The method of claim 14, including the step of 

2 repositioning a camera model responsible for 

3 taking said at least two reconstructed images. 
4 

1 17. The method of claim 16, wherein 

2 repositioning of said camera model is 

3 guided by a lookup table. 
4 

1 18 . The method of claim 17 , wherein said 

2 lookup table provides the first 

3 derivative of said translation and 

4 rotation measurements. 

5 

1 19. The method of claim 14, wherein coordinates of 

2 repositioned said camera model are transformed 

3 into coordinates for alignment of said 

4 radiation therapy beams. 
5 

1 20. The method of claim 14, wherein coordinates of 

2 repositioned said camera are transformed into 

3 coordinates of alignment of said treatment 

4 target of said patient. 
5 

1 21. A method of aligning radiation therapy beams with a 

2 treatment target of a patient, comprising the steps 

3 of : 
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4 a) acquiring a plurality of current two-dimensional 

5 images indicative of a current position of said 

6 treatment target relative to a means for 

7 generating said therapy beams; 

8 b) from a three-dimensional image comprising said 

9 treatment target, generating a corresponding 

10 plurality of reference two-dimensional images 

11 indicative of a reference position of said 

12 treatment target relative to said means for 

13 generating said therapy beams, wherein said 

14 treatment target is aligned with said means for 

15 generating said therapy beams in said reference 

16 position; 

17 c) comparing current feature vectors defined by said 

18 current two-dimensional images with corresponding 

19 reference feature vectors defined by said 

20 reference two-dimensional images; 

21 d) adjusting said current position if said current 

22 feature vectors do not match said reference 

23 feature vectors; and 

24 e) applying said therapy beams if said current 

25 feature vectors match said reference feature 

26 vectors. 

27 
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